Momentum anisotropy of the scattering rate in cuprate superconductors 
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We examine the momentum and energy dependence of the scattering rate of the high temperature 
cuprate superconductors using angle resolved photoemission spectroscopy. The scattering rate is of 
the form a + bu> around the Fermi surface for under and optimal doping. The inelastic coefficient b is 
found to be isotropic. The elastic term, a, however, is found to be highly anisotropic for under and 
optimally doped samples, with an anisotropy which correlates with that of the pseudogap. This is 
contrasted with heavily overdoped samples, which show an isotropic scattering rate and an absence 
of the pseudogap above T c . We find this to be a generic property for both single and double layer 
compounds. 

PACS numbers: 74.25.Jb, 74.72.Hs, 79.60.Bm 



I. INTRODUCTION 



There is a general consensus that understanding the 
normal state excitation spectrum is a prerequisite to 
solving the high temperature superconductivity problem. 
Angle resolved photoemission spectroscopy (ARPES) has 
played an important role in these studies because of the 
unique momentum and energy resolved information it 
provides. This includes the observation of a dramatic 
spectral lineshape change caused by the superconducting 
transition 1 ", the large momentum anisotropy of the su- 
perconducting gap consistent with d-wave symmetry 2 * 3 ., 
an anisotropic pseudogap above T<£&, the existence of 
nodal quasiparticles in the superconducting stated, and 
"strange metal" behavior above Id 7 * 8 , in the so-called 
"normal state". This state is well described by a phe- 
nomenological model called the Marginal Fermi Liquid 
(MFL) 9 , however there is as yet no consensus about its 
microscopic origin. 

A number of studies have been conducted in order to 



better understand the normal state. Of particular in- 
terest is the determination of the scattering rate, which 
governs the transport properties. Earlier on, Valla et al. 
showed for optimally doped P^S^CaC^Os+a (Bi2212) 
samples that the scattering rate measured at the chemi- 
cal potential is highly anisotropic around the Fermi sur- 
face and it varies linearly with temperature 10 . This 
anisotropy was attributed to scattering by out of plane 
impuritiesAi* 1 ^. In contrast, studies by Bogdanov et al. 13 
and Yusof et al. 14 concluded that for overdoped Bi2212 
samples, the width of the energy distribution curves 
(EDCs) (related to the scattering rate) at the chemical 
potential is isotropic around the Fermi surface. Although 
these studies of the scattering rate around the Fermi sur- 
face provide useful information, they are limited in their 
scope. They measure or infer the scattering rate (ImS) 
at a single energy (u>) and do not extract its functional 
form. Other studies that investigate the scattering rate 
as a function of energy are limited to the nodal direction 
along 13 : 14 : 15 . 

The data presented in this paper provide a comprehen- 
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sive measurement of the functional form of the scatter- 
ing rate as a function of energy and momenta around the 
Fermi surface. Consequently we are able to shed light 
on the nature of the excitations. We have conducted 
measurements of the scattering rate in under, optimally 
and overdoped samples and examine how it relates to the 
anisotropy. The key findings are: 1) In under and opti- 
mally doped samples, which show an anisotropic scatter- 
ing rate around the Fermi surface, there is a linear depen- 
dence of the scattering rate as a function of energy (for all 
points around the Fermi surface); 2) In contrast, for over- 
doped samples there is no anisotropy around the Fermi 
surface and the scattering rate is more consistent with a 
superlinear behavior with energy; 3) We also find that the 
presence of the scattering rate anisotropy in under and 
optimally doped samples and its absence in overdoped 
samples applies equally well to both single and double 
layered compounds, which suggests that bilayer splitting 
is not responsible for the anisotropy; 4) For under and 
optimally doped samples, the momentum dependence of 
the pseudogap mirrors the momentum dependence of the 
scattering rate, which implies that both may be caused 
by the same interaction; 5) The functional form of the 
scattering rate for under and optimally doped sample 
is a + buj, where the inelastic parameter, &, is isotropic 
around the Fermi surface; and finally, 6) We show that 
the bare Fermi velocity can be directly obtained from 
ARPES. 
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FIG. 1: Energy distribution curves (EDCs) along the Fermi 
surface: a) EDC data from optimally doped Bi2212 (T C =90K) 
at T=140K. b) the same for overdoped Bi2201 (T c ~ 0) at 
T=40K. The inset shows the color coded points on the Fermi 
surface where the EDCs shown in a) and b) were measured (N 
is the node and A is the antinode of the d-wave gap), c) Com- 
parison of the EDCs at the antinode for optimally and over- 
doped Bi2201 at T=40K. d) the same (at the bonding Fermi 
surface) for optimally and overdoped Bi2212 at T=100K. 



II. EXPERIMENTAL DETAILS 

Most samples employed for this work are single crys- 
tals grown using the floating zone method. They in- 
clude optimally doped Bi2212 samples (T C =90K) (used 
in an earlier study 1 ^), as well as heavily overdoped 
(T c ~ 0) Bii. 80 Pbo.38Sr 2 . iCu0 6 _5 (Bi2201) 16 . The over- 
doped Bi2212 and optimally doped Bi2Sri.6Lao.4CuOj / 
thin film samples (the latter referred to as optimally 
doped Bi2201) were grown using an RF sputtering tech- 
nique. The samples were mounted with T — M parallel 
to the photon polarization^ and cleaved in situ at pres- 
sures less than 2-10 -11 Torr. Measurements were per- 
formed at the Synchrotron Radiation Center in Madison 
Wisconsin, on the Ul undulator beamline supplying 10 12 
photons/sec, using a Scienta SES 200 electron analyzer 
with an energy resolution of 16 meV and momentum res- 
olution of 0.01 A^ 1 for a photon energy of 22 eV. 



III. SCATTERING RATES AND 
RENORMALIZED FERMI VELOCITY 

In Figs, la and b, we plot energy distribution curves 
(EDCs) along the Fermi surface in the pseudogap state 
of optimally doped Bi2212 and the normal state of highly 
overdoped Bi2201, respectively. These data reveal that in 
the optimally doped case, there is a strongly anisotropic 



pseudogap which is zero in an arc centered at the node 
of the d-wave superconducting gap (T — Y Fermi cross- 
ing), and takes its maximal value at the antinode (M — Y 
Fermi crossing). Moreover, there appears to be a strong 
anisotropy of the scattering rate, since the spectral peaks 
at the antinode are much broader than at the node. Al- 
though this has been suggested to be due to an unresolved 
energy splitting caused by bilayer mixingi^, in Fig. lc, 
we show data at the antinode for optimally doped single 
layer Bi2201, which has similar spectral characteristics 
to that of Bi2212 (shown in panel d), arguing against a 
bilayer effect as the origin of the anisotropy. We can con- 
trast the behavior shown in Fig. la with that of heavily 
overdoped Bi2201 in the normal state, where no energy 
gap is present (Fig. lb). In this case, the spectral peak 
is isotropic around the Fermi surface, indicating that the 
scattering rate is also isotropic. A similar conclusion was 
reached in recent studies of heavily overdoped Bi2212 
samples where strong bilayer splitting is presenli^^. 

To obtain more quantitative information, we analyze 
momentum distribution curves (MDCs)&. As we pointed 
out earlieji^, the MDC halfwidth (in the absence of an 
energy gap) is equal to the imaginary part of the self- 
energy at that energy, /mS(w), divided by the bare Fermi 
velocity, vpo ( n ot the renormalized one, vf). Also, the 
peak position of the MDC (as a function of binding en- 
ergy) gives the band dispersion. 
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FIG. 2: Quality of the fitting to the Momentum Distribution 
Curves (MDCs). Data for an optimally doped Bi2212 sample, 
same as in Fig. la. Each curve corresponds to a different 
energy (u>) and is measured along the k y direction at £^=0.59: 
a) lorentzian fits to the data, b) gaussian fits to the data. 



First we examine the quality of the fits to the MDCs. 
Figure 2 shows sample data from Bi2212 with lorentzian 
(panel a) and gaussian (panel b) fits. The lorentzian 
fits are a much better representation of the data at all 
energies, with a chi-squared statistic several times lower 
than for the gaussian fits. 

In Fig. 3, we plot the dispersion, Fermi velocity and 
peak widths obtained from the MDCs along selected cuts 
in momentum space, parallel to the M — Y direction. 
The first interesting point to note is that the slopes of 
the MDC dispersion curves are fairly similar in optimally 
doped Bi2212 (panel a) and vary much more in overdoped 
Bi2201 (panel b). The slope of the MDC dispersion at the 
chemical potential is equal to the renormalized Fermi ve- 
locity. We extract this from the MDC dispersion and plot 
it in panels (c) and (d) . The Fermi velocity for optimally 
doped Bi2212 (panel c) is roughly isotropic, a conclusion 
also reached in an earlier ARPES stud}^. Since in the 
nodal region the cuts along which the data were obtained 
are not perpendicular to the Fermi surface, we also show 
the velocities that were corrected by taking into account 
the angle of the cut versus the Fermi surface normal (solid 
circles). We now examine the renormalized Fermi veloc- 
ity in heavily overdoped Bi2201 (panel d). In contrast 
to optimally doped Bi2212, this is highly anisotropic and 
is much smaller at the antinode; a result consistent with 
a previous tight binding fit to normal state ARPES dis- 
persions in overdoped Bi2212*2S. The strong anisotropy 
observed in Bi2201 arises from the close proximity of 
the saddlepoint of the band dispersion at (it, 0) (where 
the velocity is zero) to the Fermi energy for this heav- 
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FIG. 3: The data obtained from MDCs (same as in Fig. fa and 
fb) for selected momentum cuts parallel to M — Y (labeled 
by the k x value along V — M, see inset to Fig. fb). (a,b) 
Band dispersion. (c,d) Fermi velocity - raw data are shown 
as empty circles, the data corrected by the angle of the cut 
are shown as dots. (e,f) peak widths, a), c), e) optimally 
doped Bi22f2 at T=I40K and b), d), f) overdoped Bi220f at 
T=40K. 



ily overdoped sampled. Stated another way, the small 
anisotropy in the optimally doped case implies that the 
saddlepoint in the band structure at (7r, 0) is significantly 
far away from the Fermi energy. 

This effect is illustrated in Fig. 4a. Here we plot the 
band dispersion obtained from a tight binding fit where 
the saddlepoint has been assumed to lie far below (150 
meV) the chemical potential - the same energy as at (7r, 0) 
in optimally doped Bi2212. The band dispersion from the 
tight binding fit is plotted along cuts taken in the same 
geometry as the experimental data (Fig. 3). In panel 
(b) we show the Fermi velocity calculated from the slope 
of the band dispersion. This velocity bears a close re- 
semblance to the real data for optimally doped Bi2212 
in the sense that is roughly isotropic (Fig. 3c). We now 
address panels (c) and (d) , which show the tight binding 
dispersion and Fermi velocity, respectively, in the case 
where the saddlepoint at (ir, 0) lies near (10 meV below) 
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FIG. 4: The band dispersion (a,c) and renormalized Fermi 
velocity (b,d) obtained via a tight binding model where the 
saddlepoint at (-7T, 0) lies at (a,b) 150 meV below the chemical 
potential (equivalent to optimally doped Bi2212) and (c,d) 10 
meV below (equivalent to overdoped Bi2201). 



the chemical potential - similar to the situation in over- 
doped Bi220li&. The effect of the saddlepoint is readily 
apparent as the velocity decreases quite rapidly near the 
(tt, 0) point - exactly the same as in the experimental 
data (Fig. 3d). In summary, our results for the Fermi ve- 
locity are consistent with tight binding calculations and 
suggest that renormalization effects in the normal state 
are momentum independent. This strongly supports our 
conclusion (see later) that the b term of the self energy 
in under and optimally doped Bi2212 is also momentum 
independent (isotropic). 



IV. FUNCTIONAL FORM OF THE 
SELF-ENERGY 

The energy dependence of the MDC peak widths for 
the various momentum cuts is shown in Figs. 3e and 3f. 
Consider first the optimally doped Bi2212 data (3e). To 
a good approximation, the result for any particular cut 
can be fit to the form om + bu^^- This is analogous 
to the am + bjtfT form indicated for the temperature 
dependence of MDC widths previously reported by Valla 
et aliS.. 

In Fig. 5a, we show the momentum dependence of the 
a,M and 6m terms extracted from the MDC HWHM for 
optimally doped Bi2212 in Fig. 3e. To complement these 
results, we have also fit the EDCs along the Fermi sur- 
face for this sample using a model self-energy, E. We have 
tested both quadratic and linear energy dependences for 
JmE and found that only the latter gives an adequate 
description of the data. i?eS is determined by Kramers- 
Kronig transformation of JmS = a + bui, assuming the 
latter saturates at a constant value beyond a cutoff en- 
ergy of 0.5 eV. The background signal was determined us- 
ing methods described previously 22 and added to the cal- 
culated curves. The experimental energy resolution was 
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FIG. 5: The pseudogap and the elastic ("a") and inelastic 
("6") portions of the scattering rate (ImT,—a + bui) around 
the Fermi surface for optimally doped Bi2212. a) momentum 
dependence of the aj\/ and &m terms obtained by fitting MDC 
HWHM data from Fig. 3e. b) momentum dependence of the 
a and b terms obtained by fitting EDCs from Fig. la. c) 
bare velocity obtained by dividing the b and 6m coefficients - 
raw data are shown as empty circles, the data corrected by the 
angle of the cut are shown as dots, d) Position of the midpoint 
of the leading edge of the EDC around the Fermi surface for 
optimally doped Bi2212 obtained from Fig. la. This is an 
approximate measure of the pseudogap - the actual value of 
the pseudogap is equal to about twice the midpoint shift—. 
(k x labels the momentum cut as in Fig. 3, with k x = 0.4 
corresponding to the node and k x = 1.0 to the antinode). 



taken into account by convoluting the calculated curves 
with the appropriate gaussian function. The experimen- 
tal momentum resolution was taken into account by sum- 
mation of the spectra over the analyzer momentum win- 
dow. The quality of the fits is illustrated in Fig. 6 along 
with plots of the relevant self energies S. To improve the 
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determination of the b coefficient we have performed fits 
to EDCs peaked at high binding energy, where b has the 
biggest impact on the lineshape. The EDC fits performed 
close to the Fermi momentum allow a precise determina- 
tion of the a coefficient. In Fig. 5b, we show the values of 
the a and b coefficients obtained from these EDC fits 23 . 
Note the similarity of the results to those in panel a) de- 
spite the quite different methodologies used. Hence, we 
can have confidence in the validity of our results for the 
a and b parameters. 



V. EXTRACTION OF THE BARE FERMI 
VELOCITY FROM ARPES 

As an interesting aside, we can exploit the above to 
estimate in a very simple way the bare velocity. We have 
previously showni^ that the width of the MDC peak is 
given by Wm = £"(w)/«fo- Since now Wm — im+^mw 
and ImH=a + bu>, it follows that vfo = b/bM, therefore, 
we can extract the bare velocity around the Fermi surface 
directly from our data as shown in Fig. 5o~. We note 
that the bare velocity obtained in this simple way is con- 
sistent with band theory predictions and is in agreement 
with alternative method of extracting bare velocity 2 ^. 



VI. CONCLUSIONS 

Returning to Fig. 5a and 5b, the first conclusion to 
draw from the data is that the b (6m) term is isotropic 
in both plots. At first sight, it would appear that the 
Bi2201 case is different, since the &m term in that case 
(slope of the curves in Fig. 3f ) appears to increase as the 
antinode is approached. But once the velocity is divided 
out (Fig. 3d)2£, we find in this case as well that the 6 term 
for JmE is isotropic, which is consistent with the isotropy 
of the EDC lineshapes shown in Fig. lb. The isotropy of 
b provides strong support of the original marginal Fermi 
liquid conjectural, where isotropic (i.e., local) behavior 
is required to guarantee uj/T scaling. 

The om term (obtained from the zero intercept in 
Fig. 3e) is highly anisotropic. This is consistent with 
the strong anisotropy of the EDC lineshapes shown in 
Fig. la. The anisotropy of the a term in optimally doped 
samples has been attributed to off planar irnpuritiesii*^. 
On the other hand, we note the remarkable similarity 
between the anisotropy of this term and that of the 
pseudogap 2 ^ (Fig. 5d). This indicates to us that the 
anisotropy is probably not due to impurity scattering, 
but rather is related to the same interaction that gives 
rise to the pseudogap. This is consistent with the obser- 
vation of isotropic lineshapes for more heavily overdoped 
samples of Bi2212 where no pseudogap is present. How- 
ever, other possibilities could also be considered, such as 
the cold spots model of Ioffe and Millis, where a highly 
anisotropic scattering rate is conjectured due to scatter- 
ing from d-wave pairing fluctuations 29,30 . 
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FIG. 6: Fits to energy distribution curves using an a+bcj 
model for the self energy, a) real and imaginary parts of the 
self energy used to fit the EDC data in the nodal region, b) 
real and imaginary parts of the self energy used to fit the 
EDC data in the antinodal region, c) EDC data and fits in 
the nodal region. The upper curves are close to the band min- 
imum (k y =0), while the lower curves are close to the Fermi 
momentum, d) EDC data and fits in the antinodal region. 
The upper curves are close to the band minimum (k y =0), 
while the lower curves are close to the Fermi momentum. 



Finally, we note that an anisotropy in the "zero inter- 
cept" of the MDC width is also evident in the heavily 
overdoped Bi2201 sample (Fig. 3f). But, in this case, it 
can be accounted for by the anisotropy in the bare Fermi 
velocity, which decreases towards the antinode. Thus the 
a term in JmE is isotropic (again, this is consistent with 
the isotropy of the EDC lineshapes shown in Fig. lb). 

In conclusion, we find that the normal state scattering 
rate in under and optimally doped cuprates can be ap- 
proximated by the form a + bu). The inelastic b term is 
found to be isotropic, which is a necessary ingredient in 
the marginal Fermi liquid conjecture 9 -. In contrast, the a 
term is found to be anisotropic for under and optimally 
doped samples, with the anisotropy linked to that of the 
pseudogap. 
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